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ABSTRACT 



We argue that any violent galactic winds following early epoch of star bursts would signif- 
icantly weaken the potentials of galaxies, and leave lasting signatures such as a lowered dark 
halo density and preferentially radial/escaping orbits for halo tracers such as globular clusters. 
A galaxy is disintegrated if more than half of its dynamical mass is blown off. The presence of 
dense halos in galaxies and the absence of intergalactic/escaping globulars should imply an upper 
^ ■ limit on the amount of baryons lost in galactic winds around 4% of the total mass of the galaxy. 

I This translates to limits on the baryons participating the early star bursts and baryons locked 

' in stellar remnents, such as white dwarfs. The amount of halo white dwarfs claimed in recent 

I proper motion searches and microlensing observations in the Galactic halo are too high to be 

' consistent with our dynamical upper limits. Similar arguments also imply upper limits for the 

ps) I amount of neutron stars and stellar black holes, in galaxy halos. Nevertheless, a milder outflow 

' is desirable, especially in dwarf galaxies, both for lowering their cold dark matter central density 

and for injecting metals to the intergalactic medium. 

I ' Subject headings: Galaxy: halo - Galaxy: kinematics and dynamics - galaxies: interactions - dark matter 

O ' - galaxies: clusters 

^ 1. Introduction 



■ How strong is the energetic feedback from star formation in galaxies? What are the effects of massloss 

, from a galaxy? These are important questions in understanding galaxy formation. Baryonic gas cools and 
' forms stars in potentials of the dark matter halos. Some massive stars explode as supernovae and the stellar 
envelope is ejected with high speed. This feedback not only injects and mixs metals in galaxies, but also 
drives galactic winds to pollute intergalactic medium. A sudden loss of baryonic gas by winds after an 
extremely powerful burst of star formation also means a sudden weakening of the gravitational potential of 
the galaxy. 

The weakening of the potential has several interesting consequences. For example, the galaxy halo will 
relax to a looser distribution after a mild loss of its mass through winds. This reduction of halo density 
is interesting for overcoming the dense cusp in Cold Dark Matter halos (e.g. Gnedin & Zhao 2001). Stars 
could also escape a lowered potential well of a galaxy, and appear as intergalactic stars. 

There are also limits on the amount of massloss. A galaxy becomes unbound if more than half of its mass 
is lost (Hills 1980). Real galaxies appear to be tightly bound. The boundness of present day galaxies should 
set some limits on the violentness of massloss and star formation in the past. 

It is now well-established that galactic outflows are ubiquitous in actively star forming galaxies, both in 
the local universe and at high redshift (see a review by Heckman 2001). The amount of gas lost in the 
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energetic outflow is largely comparable to the amount of stars formed. Such a wind can be efficient in 
polluting the hot gas in the intergalactic medium (IGM). Some amount of energetic wind or outflow seems 
mandatory to account for the presence of metals in the QSO absorption systems at high redshift, and metals 
in the hot gas in galaxy clusters (e.g., Bookbinder et al. 1980). 

In the Milky Way there are very few metal-poor stars with [Fe/H] < —4 in the halo and [Fe/H] < —0.6 
in the disk. There are also metals in some distant high-velocity clouds in the Local Group. These again 
suggest an early phase of metal production in a burst of mainly massive stars. 

Several recent claims of detection of a substantial amount of white dwarfs (WDs) in the Galactic halo and 
distant galaxies by direct proper motion searches or indirect microlensing observations suggest a first burst 
of high- mass stars narrowly peaked around (2 — 3) Mq in galaxy halos. To leave behind every half-solar-mass 
WD, a few times more of the remnant mass is returned to the interstellar medium via the planetary phase. 
Those in binaries die after ~ 1 Gyr via the highly energetic Type la SN, which can easily power a hot and 
fast wind and drive all stellar ejecta out of the potential well of a galaxy. 

Motivated by these considerations, wc model the dynamical effect of the post-star- formation galactic wind 
on the density distribution of the dark particle halo. We examine the effect of a fast wind at a redshift z ~ 1 
when the Milky- Way-sized halos are already largely assembled, and energetic SN la are observed. In general, 
a galaxy will expand, stellar orbits become highly radial in response to a rapid weakening of the potential well 
due to the mass gone with the wind. In particular we set limits on the amount of high-mass star formation 
in galaxies in the past, and the amount of stellar remnents, e.g., white dwarfs, in present day galaxies. 

Simple models for the wind and expansion are developed in §2 and the Appendix, including the effects 
on globular clusters. Applications to the halo white dwarf problem are made in §3. We discuss our results 
in §4, in particular, the strength and starting time of the wind, and conclude in §5. 



2. Models 



2.1. The post- wind halo expansion and escaping stars 

Let us define e as the fraction of the initial galaxy total mass that is lost in the wind. It is then given by 

where Mq is the present day total mass bound to the galaxy, and My, is the amount of mass lost in the 
wind. Define the constant eoo to be the critical amount of massloss for the system to expands to infinity or 
unbound, then for a wind of a range of speed we have 

[fast wind] 0.5 < Cqo < 1 [slow wind], (2) 

where the boundaries are specified such that a system unbounds if lossing 50% of its mass in a fast wind or 
100% of its mass in a slow wind. 

The loss of gas reduces the gravitational mass of the system, making the potential shallower. The system 
will relax to a new equilibrium by an overall expansion. The orbits of tracer particles in the galaxy potential 
will also expand. The expansion factors depend on the speed of the massloss, and are derived in Appendix 
A.l and A. 2. The weakening of the potential well also means some stellar particles could escape. The process 
is different from the massloss of the gas, which is pushed out by the pressure gradient. Here stars at the high 
energy tails could gain sufficient potential energy and exceed the bounding potential of the system. This 
fraction is estimated in Appendix A. 3. 

In general a system of initial characteristic size ri,e expands, after lossing e fraction of the total mass, to 
a size r/,e with 

X = - = 1 + ^22^. (3) 

Ti Coo C 

Note X — > oo when e e^o- For a slow wind defined by = 1 we have X = (1 — e)~^. 
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The weakening of the potential due to massloss and overall expansion has also the effect of increasing 
the apocenter of a tracer orbit. A tracer particle initially on circular orbit of size rj = Oj will expand to a 
generally eccentric orbit with an apocenter r = a, where the expansion factor 

Coo 

so in a fast wind ^ the expansion in apocenter is somewhat stronger than for the system as a whole, i.e.. 
A = — > — — X . This is because the pericenters are determined by the angular momentum barrier, which 
does not change very much, so the apocenter must increases faster than the semi-axis. As a result, orbits 
become preferentially radial. Note also the size of a tracer orbit doubles, i.e., ^4 = 2 for e = Cdbi = O.Scoo- 

Another effect of weakening of the potential is that some stars or globular clusters can escape the potential. 
Let e* be the fraction of escaping stars among all stars produced in the galaxy, then (cf. Appendix A. 3) 

2 

(5) 



(1 - eoo)e 



(1 - e)eo 



Note that no stars escape in a slow wind with Cqo ^ 1, and all stars escape in a wind of maximum strength 
e ^ Coo- If the massloss is small, then e* in a fast wind with Cqo = 0.5. 

2.2. Absence of intergalactic globulars and the maximum massloss 

We are primarily interested in fast winds with e^c = 0.5. We shall also assume that the expansion factors 
X and A and the massloss fractions e and are all everywhere the same. These arc illustrated in Fig.l 
as functions of gas massloss fraction e. Since galactic winds are generally expected to be fast, the above 
calculations motivate us to rate the degree of expansion caused by a fast wind into three categorical classes, 

e ~ 0.04 mild, (6) 
0.15 severe, (7) 
~ 0.5 unbound. (8) 

To justify these categories, note when a galaxy losses e = 15% — 25% of its mass in a fast wind, the galax;y 
expands in volume by a factor X^ = 1.8 — 3.4, and the apocenters of the globular cluster tracers could 
increase by a factor of A ^ 1.4 — 2, and a fraction = 3% — 10% of the galactic stars or globular clusters 
should escape. For a bright galaxy like the Milky Way or the Andromeda with of order 150-450 globulars, 
this would predict of order 5-45 "intergalactic" globulars. 

Fig. 2 shows the most recent compilation of orbital parameters of the Milky Way globular clusters (Dinescu 
et al. 1999 and references therein). The distribution appears similar to that of an isotropic models, which 
predict a median apo-to-peri ratio about 4 (van den Bosch et al. 1999). Globular clusters with apo-to- 
peri ratio > 10 are rare, and these outliers appear as frequently for inner disk globulars as for outer halo 
globulars; the outer globulars should experience much less dynamical friction and tidal destruction than the 
inner globulars. There are very few Galactic globulars beyond 20 kpc of the Galactic center. The total 
number of globular clusters in the Local Group is estimated to be 

A^MW + NAnd 150 + 450 = 600 (9) 

(e.g., Barmby & Huchra 2002), and the numbers might have been much higher earlier on considering the 
destruction of inner globulars once the disk and bulge form. The fact that there are no direct evidence for 
free floating clusters in the Local Group suggests at face value 

^For a slow wind we have eoo = 1, hence A = X = 1/(1 — e) (cf. eqns. 3 and 4), consistent with the expectation that the 
eccentricity of the orbit is constant for adiabatic changes of the galaxy potential (Lynden-Bell 1963) 
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And as far as wc know there are no direct evidence for intergalactic globulars around other galaxies, implying 
perhaps that they are very few. Together these imply that the massloss in galaxies in the past is limited to 
a fraction 

e < 0.04. (11) 

We conclude the lack of escaping old halo globular clusters in galaxies does not support a rapid expansion 
of galaxy halos. 

3. Application to the halo white dwarfs 

3.1. The observational and theoretical case for halo white dwarfs 

Recent searches for baryonic dark matter found between (2-5)% to (8-50)% of the dark halo density to 
be in the form of half-solar-mass WDs. The lower estimates come from direct detection of faint blue high 
proper motion objects in the solar neighbourhood (Oppenheimer et al. 2001, Ibata et al. 2000), and the 
upper estimates are indirectly inferred from observed microlenses towards the LMC (Alcock et al. 1997, 
2000, Lasserre et al. 2000). While these observations confirms the notion that the bulk of the halo is still 
in the form of undetected non-baryonic weakly interacting particles (WIMPs), the existence of a significant 
halo WD population is still very contentious (see a recent review by Richer 2001). In fact the infered 
amount of WDs is surprisingly big, and is difficult to fit in current theories of star formation and galaxy 
evolution. For example, the mean baryons to WIMPs split in the universe is only around 15 : 85 for current 
favored cosmological parameters. Nevertheless, theorists have long been prepared to include this odd piece 
in the cosmic puzzle (Carr 1994). The most natural way to account for these WDs and their extragalactic 
counterparts, albeit not without difficulties, is that they are the relics of the first burst of high-mass stars in 
all galaxies. 

The detailed properties of the burst of WD progenitors must conform with a number of observational and 
theoretical constraints. The progenitors are likely formed about 10 Gyrs ago, earlier than the disk, so that 
the hydrogen-rich WDs have enough time to cool to a temperature of 3000-8000 K after their progenitors 
evolve off the main sequence in 1-2 Gyrs (Hansen 1999); the helium- rich WDs can cool faster, and are perhaps 
as numerous as the hydrogen-rich WDs. The progenitors must have a narrow top-heavy IMF within the 
range of 1 — 8 Mq, perhaps a log- normal distribution peaked around 2 — 3 Mq, so to avoid overpredicting the 
counts of faint main sequence stars in the halo (Chabrier, Segretain & Mera 1996 and references therein). 
This burst of WD progenitors must have a total mass at least comparable to the total mass of galactic stars 
in the disk today, even using the most conservative lower estimate of Oppenheimer et al. (2001). Note we 
treat the formation of WD progenitor as a single burst, although some claim a spread in ages among the 
sample of white dwarfs from Oppenheimer et al. (Hansen 2001). 

The bulk of the mass of a WD progenitor is in the envelope, which is injected back into the interstellar 
medium at the end of main sequence lifetime through the planetary nebulae phase. The metals synthesized 
in these progenitors are also released, which could severely distort the abundance pattern of observed halo 
and disk stars, especially in deuterium, helium, carbon and nitrogen, if the metals were kept inside the 
galaxy. This is saved by a galactic wind about 1 Gyr later powered by Type la supernovae (SN la ) evolved 
from binary progenitors (Fields, Mathews, Schramm 1997, Fields, Freese, Graff 2000 and references therein). 

3.2. A dynamical limit on white dwarfs 

Let us estimate the effect of a wind driven by such a burst of WD progenitors. Let Mq be the present day 
galaxy mass, /wdMq be the mass locked in WD remnants. For a progenitor mass of mpc to evolve into a 
WD of mass mwD, the amount of mass lost is mpo — mwD, so the total number of WDs (as well as the 
WD progenitors) N\yd and the mass lost in the wind is given by 

M^ = Mi-Mo = {mpG - mwo) Nwd, (12) 
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and 

N^n = (13) 
mwD 

where M, = Mq + is the pre-wind original total mass of the galaxy. The fraction of the initial galaxy 
total mass that is lost in the wind, e, is then given by 

( mpg—mwn \ t 
^^r. )f^^ ^^^^ 



mwD 



WD 



Fig. 3 shows the boundaries defined by eq. (14) and eqs. (6,7,8) in the parameter plane of mpc/mwD vs. 
fwD- For a log- normal distribution of mpc within 1 Mq — 8 Mq, we expect 

(mpc) - 4 (mwD) - (2.5 - 4.5) Mq, (15) 

insensitive to tlw^ details of the IMF. So, allowing only a mild wind and expansion. Fig. 3 predicts that the 
maximum mix of halo WDs with halo WIMPs is given by 

fwD . 2% 



1-fwD- 98% 
in all tightly bound galaxies. 

4. Discussions 

4.1. Did the Galaxy experience a severe expansion? 

Apart from the globulars, another constraint for the amount of evolution of the Galactic potential is that 
the present-day dynamical mass of the Milky Way out to a shell at 50 kpc (the distance of the LMC) is about 

0. 5 X 10^^ Mq (Wilkinson & Evans 1999), implying an average density of ~ 10~^ MqPC"^ inside this shell. 
One can also estimate the density of an initial dark halo of the A4ilky Way, averaged over the same radius. 
Assume the inner 50 kpc of the Galaxy is already virialized at z = 2, then it should have an overdensity of 
at least some 200 times the mean density of the universe then, insensitive to cosmology. More precisely a 
collapsed halo has a mean density 

Pvir(^) ~ lM^t^£M(i + ,)3 ^ 0.8 X 10-3 MePC-3, (17) 

where Ho = 75km/s/Mpc, z = 2 and flz ^ 0.9 for flo = 0.3 and fl^ = 0.7 cosmology (e.g., Bullock et al. 
2001). So the density of the dark halo at z = 2 agrees with the present day halo within 50 kpc very well. 
Another way to estimate this is with the Navarro, Prenk & White (1996) or Moore et al. (1998) models, 
assuming a virial radius of the Milky Way halo (200 kpc-400 kpc), and the concentration parameter of 6-12. 
These models can predict the mean density within 50 kpc, and this turns out to be (0.5 — 2) x 10"^ Mqpc-^, 

1. e., within a factor of two of current density. In comparison if a wind doubles the size the galaxy, the mean 
density enclosed inside a shell drops by a factor of = 8. This implies that any dynamical evolution must 
have been very mild after the formation of the outer halo. 

4.2. When did the galactic wind start? 

The oldest halo globular clusters are as old as about 12 Gyrs (Salaris et al. 1997), and the coolest local 
disk WDs are 9-12 Gyrs (Knox et al. 1999, Oswalt et al. 1996, Hernanz et al. 1994). It is plausible 
to assume that 9 — 12 Gyrs ago is roughly the epoch when baryons collapse in the potential of the dark 
particle halo of the Milky Way, and start forming the first stars everywhere, in the local disk, in globular 
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clusters and in the halo; the IMF is top-heavy in the outer halo and normal Salpeter elsewhere in the disk 
and in globular clusters. This epoch corresponds to a redshift of z ~ 1 — 3 in a 14-Gyr-old universe with 
1 — flm = fl\ = 0.66 and Ho = 66km/s/Mpc. The Madau diagram, i.e., the extragalactic star formation 
rate in this high redshift window is still not yet well constrained. There is a delay between the onset of the 
halo star formation and explosion of SN la , which is the time for the secondary of a binary to evolve and 
overflow its Roche lobe, and bring the WD of the primary into a critical mass. Take a plausible mass of 
^ (4 — 8) M0 for the primary, and ^ (1 — 2) M0 for the secondary, the delay is ^ 1 — 8 Gyr. This predicts 
the appearance of the first halo SN la ~ 8 — 11 Gyr ago, which seems to be consistent with current detection 
of some SN la at redshift z = 1 — 2. In short, we estimate that the galactic wind powered by the first SN 
la in the halo starts at z ^ 1, or when the universe was about half its present age. Since major mergers 
are rare in Milky Way sized halos after z = 2 in the fairly successful ACDM cosmology, it appears safe to 
assume that CDM halos as big as the Milky Way halo are well assembled through hierarchical merging and 
infall by the time of the wind. 

4.3. Are the winds fast or slow? 

A slow wind means nearly adiabatic expansion, which happens if the outflow velocity is much less than 

the circular velocity of the galaxy. Observationally prominent outflows are seen in active star burst galaxies 
such as M82 at low redshift. Generally speaking the outflow rates are comparable to the rates of star 
formation in terms of mass and energy output, and these superwinds driven by star bursts have velocities of 
100 — 1000 km s"""^ (Heckman 2001). A wind or jet must also carry sufficient momentum in order to deliver 
and mix the metals observed in the intergalactic medium. At high redshift winds with Doppler offset velocity 
of 200-1000 kms~^ have been seen in Lyman-break galaxies with active star formation (e.g., Pettini et al. 
2002), suggesting for fast outflow driven perhaps by UV- photons and SNe from the star burst. It would be 
oversimplifying to associate these observed star bursts directly with the speculated halo WD progenitors, 
but it is perhaps safe to argue that a hypothetic burst with a top-heavy IMF should lead to outflows at least 
as violent as these observed ones. 

We could also estimate the energy budget of winds theoretically, as in Fields, Freese & Graff (2000). 
Ejecta from a SN la has a velocity v^j > ^lO^iergs/8 M© ~ 2500 km s~^, where 10^^ erg is the kinetic 
energy per SN la , and 8 M© is the upper limit for the mass of the exploded star. So the ejecta could escape 
any giant galaxy instantaneously by its own kinetic energy. It would require fine tuning of the star formation 
efficiency and/or the binary fraction to keep the outflow velocity low, at least for a top-heavy IMF.^ As an 
illustration, consider a very inefflcient star burst early on with every 1000 Mq of gas of a galaxy forming 10 
single stars of 6 M© plus one binary of masses 4 M© : 2 M©. Assume this top-heavy IMF leaves behind a 
remnent of a 1 M© WD with another 5 M© of stellar material ejected back to the interstellar medium with 
a velocity of SOkms""'^ for a single progenitor, and 3000 km s""'^ for a binary progenitor. Most part of the 
ejecta might escape the galaxy with a high speed through "holes" , however, if only 25% of the ejecta could 
mix and thermalize with another 934 M© of unprocessed gas, then the mean outflow velocity is still above 
3000 km s"^ ^25% x 5 Mq/989 M© ^ 100kms"\ i.e., the outflow velocity is stiU far from adiabatic with a 
speed comparable to the velocity dispersion (sound speed) for the Milky Way-sized galaxy, and exceeding the 
escape velocity of dwarf galaxies. So an adiabatic expansion in the Milky Way would require "tuning" the 
binary formation efficiency <C 10%, the star formation efficiency <C 6.6% and the energy deposit efficiency 
<C 25%. Without this fine-tuning, we expect theoretically that the energy from the SNe is too high for a 
slow galactic wind. ^ 

In short, we conclude from theoretical and observational evidences that the galactic winds have speeds 

^Using a bottom-heavy IMF can suppress (binary) stars above 1 Mq and the SN la , but it is not interesting for producing old 

halo WDs in the first place. The ejecta of low-mass stars drives a wind of velocity of ~ 30kms~^, fast for dwarf galaxies. 
*It is a common misconception that pressure-driven winds have necessarily a sound velocity very close to the typical velocity 

dispersion of the dynamical system. Calculations of winds in globular clusters show that the sound velocity can exceed the 
typical velocity dispersion and escape velocity of the dynamical system by a factor 1-10 if the energy injection rate is high (see, 
e.g.. Fig. 2 and Fig. 3a of Faulkner & Freeman 1977). The gas flow velocity can also exceed the sound velocity in some regions. 
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somewhere between our slow wind and our fast wind limits, with a tendency to be fast winds. For simplicity 
we shall adopt the fast wind limit in this paper, which could somewhat overestimate the effects of the wind. 

4.4. Comparison with independent observations 

The most direct estimates for the WD fraction fwD come from direct searches by Ibata et al. (2000), and 
by Oppenheimer et al. (2001); the latter claims a firm lower limit from their large sample of high proper 
motion halo WDs, about 38 with spectroscopy confirmation and some show hydrogen lines, indicative of 
a mix of hot and cold WDs with temperature in the range 3000 — 8000 K. Indirect measurements come 
from accounting for 10 — 20 microlensing events towards the LMC from the MACHO group (Alcock et al. 
1997, 2000) and from the EROS group (Lasserre et al. 2000); the latter claims only an upper limit on the 
basis of a null detection of a conclusive halo microlens. Clearly our theoretical model is consistent with the 
EROS upper limit, but inconsistent with the best values from the MACHO group in 1997 and 2000. It is 
possible that most of the microlensing towards the LMC are dominated by stellar microlenses in the LMC 
(Sahu 1994, Zhao & Evans 2000) with only a minor contribution from the halo WDs. Our theoretical result, 
eq. (16), which is effectively an upper limit on any WDs in the halo, is consistent with the most recent lower 
limit by Oppenheimer et al. (2001). 

Our dynamical constraint is the tightest for the high mass progenitors. One would hope that the progenitor 
mass could be constrained by the temperature of observed WDs together with a cooling model and an 
estimation of epoch of formation. For example, the cosmic clock in Fig. 3 indicates a 10-Gyr-old burst. Their 
WD remnants are now hotter than 4000 K if nipG > 2.2 Mq, and cooler if m,pa < 2.2 Mq. Unfortunately 
progenitors of 2 Mq and 4 Mq formed 9—12 Gyrs ago both evolve to WDs in < 1 Gyr, which then cool 
to nearly the same temperatures after another ~ 10 Gyrs (cf. Fig. 4), so we cannot differentiate them. We 
could, however, separate the WDs of 1 Mq progenitors from WDs of 4 — 8 Mq progenitors Hinc;e the former 
will have a temperature (4000-10000 K) much higher than the latter (< 4000 K). The temperatures of the 
observed halo blue WDs are spread out in the range 3000 - 8000 K (Oppenheimer et al. 2001). This means 
there must have been some 1 Mq progenitors and some 8 Mq progenitors, consistent with a log-normal IMF 
peaked around 2 — 3 Mq (Chabrier et al. 1996). 

Fig. 5 shows the dynamical limits on cool and hot WDs from a 10-Gyr-old population. Note stronger limits 
on cooler WDs because these have more massive progenitors with higher mpc/'mwD ratio. For a typical 
temperature of 4000 K, we find fwD < 2%, consistent with our earlier dynamical limit (eq. 16). 

4.5. The minimal dynamical feedback 

The dynamical feedback is generally strong. One way to minimize it is to form the WD progenitors in 
a thick disk instead of in the halo. It has been argued that some of the high proper motion WDs in the 
Oppenheimer et al. (2001) sample have thick disk kinematics (Reid, Sahu fc Hawley 2001, Koopmans & 
Blanford 2001). Assuming all WDs and their progenitors are in a thick disk of a scale height 1 kpc, we can 
convert the WD number density of 2 x 10~'*pc^"^ of Oppenheimer et al. (2001) into a local column density 
of progenitors of 0.4mpG Mgpc"^. In comparison, the local column density of the Galactic thin disk is 
^ (30-40) Mqpc~^ within \Z\ < 1 kpc, and the dark halo column density is ~ (100-200) Mqpc"^ within 
a vertical column \Z\ < oo depending on the density profile of the halo. We estimate ~ (0.6 — 1.2)% of the 
dynamical mass in the local column is in WD progenitors if mpc ~ 3 Mq. Take e = 2% as the characteristic 
fraction of the dynamical mass for the inner halo lost in a wind, then the inner halo would expand by a 
factor X ^ 1.02 for an instantaneous wind. It is interesting to compare this with the dynamical effect of 
the formation of a galactic disk collapsing out of the spherical halo. The thin disk makes up ~ 15% of the 
dynamical mass in the local vertical column with \Z\ < oo. If a galactic thin disk with 15% of the mass of 
the halo form adiabaticly, the halo contracts by a factor about 0.85. So the feedback from a thick disk WD 
population is dynamically unimportant compared to the formation of the galactic thin disk. 
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4.6. Formation of baryonic disks and mild wind 

How docs the amount of mass lost in a mild wind compare with the total amount of baryons in a galaxy? 
Take a mean density of baryons in the universe Qb = 0.05(-ffo/66)~^, and the fairly successful ACDM 
cosmological models with 1 — = 0,^ = 0.66 and -ffo = 66km/s/Mpc, we expect a typical baryonic fraction 
in galaxies 

^ ~ 15%. (18) 

So we can express the wind massloss in terms of the total baryons in a galaxy. A mild wind (cf. eq. 6) looses 
25% of the total baryons of a galaxy 

~ 25%M6, Mb = (19) 

A mild wind is also relevant for star formations in galaxy disks with a normal Salpeter type initial mass 
function with very few massive stars. Here the baryons of a galaxy cool and collapse adiabatically to a disk. 
When star formation ends, the envelope of evolved stars are ejected out with a mild galaxy wind; the amount 
of stellar ejecta ^ 25% of the initial stellar mass for a bottom-heavy IMF, so the wind is mild. We will 
not digress into this, except remarking two counter-acting effects of disk formation on the density profile of 
galaxy halo: the collapse of gas into disk stars leads to contraction of the halo, and the later massloss from 
winds leads to expansion of the halo. 

4.7. Disintegration of galaxies by winds and intergalactic stcirs and globulars 

What are the chances and consequences for a heavy wind to disintegrate some galaxies in the past? 
Systematic destruction of dwarf galaxies would be desirable because of the notorious overprediction of dense 
dwarf satellites in generic CDM models. Another consequence of disruptions is that the stars and their 
remnants, and perhaps globular clusters can disperse with the unbound dark halo into the intergalactic 
space and are later recycled to into the dark halo of the second generation massive galaxies, such as the 
Milky Way. Let /& be the baryonic fraction in a galaxy, we expect in ACDM models a typical baryonic 
fraction in galaxies (/b) ~ rib/rim = 15%. To disrupt a galaxy, we need to have about half or more of a 
galaxy to be in baryons, i.e., /b > e > 0.5, where e is the wind-lost baryons as a fraction of the mass of 
a galaxy. So a wind-disrupted galaxy needs to be very baryon-rich, three times that of the universal value 
(/b) . So the chances for disruption are likely low, unless the universal baryonic fraction is much higher than 
15%. 

There are no strong evidences for intergalactic stars or globular clusters although some have been claimed 
in galaxy clusters. Both M87 and NGC1399 have excess star lights and numerous globulars around them, 
and the outer globulars have velocities which are more characteristic being bound by the galaxy cluster 
rather than the central galaxy (e.g., Kissler-Pratig et al. 1999). Theuns & Warren (1997) also discover 
10 PNe in the Fornax cluster but unassociated with any bright galaxies. One complication of interpreting 
galaxy cluster data is the projection effect: the PNe or globulars might be "markers" of halos of low surface 
brightness galaxies in the galaxy cluster in the line of sight. Even the bona fide free-floaters are likely striped 
from spiral galaxies due to frequent merging in a galaxy cluster (Moore et al. 1999), which is very different 
mechanism from the wind in an isolated galaxy. In any case the cD galaxies are ambiguous systems for 
constraining galaxy winds compared to the Local Group galaxies. 

In principle some amount of free-floating globulars cannot be ruled out in the Local Group. When placed 
at 1 Mpc, the edge of the Local Group, a 1O'''L0 globular would appear as bright as a K-dwarf star at 1 kpc 
in the halo, but has a half-light radius of about 5 pc ~ 1". Much of the sky has not yet been covered by 
detailed searches except for perhaps a few thousand square degrees: along the Magellanic Stream (10° x 100°) 
and the Sagittarius Stream (10° x 360°) and around Andromeda (10° x 10°). High quality data over 1500 
square degrees has been released by the on-going SLOAN survey. Assume 10% of the Local Group is well 
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surveyed and no new globulars are found, then < 10 globulars could be by chance at the edge of the Local 
Group in unsurveyed area. 

5. Conclusions 

In short, the absence of obvious signs of severe expansion of the Milky Way and the absence of escaping 
globulars in galaxies in general suggest that star formation and feedback in galaxies are mild; galactic winds 
carry less than 4% of the total gravitational mass, or 25% of the baryonic mass in galaxies. This restricts 
the amount of remnents in galaxy halos from early star formation, and cool halo WDs can make up 2% of 
the total mass in present day galaxies. This is consistent with only the lower end of the halo WD fraction 
of 2% < fwD < 50% from direct and indirect detections of WDs; it is consistent with Oppenheimer et al. 
(2001) detection, but inconsistent with the higher values from Alcock et al. (1997, 2000) and Ibata et al. 
(2000). Our dynamical limit is also tighter for very cool WDs than for hot WDs. There is very limitted 
room to fit in a large extra population of very cool WD population < 4000 K, something that future deeper 
surveys should take into account. 

Our model implies stringent limits on the amount of black hole (BH) or neutron star (NS) formation as 
well. Assuming these are remnents of massive stars formed in galaxy halos, then typically neutron stars 
or stellar black holes are formed from very massive (8 M0 — 1000 Mq) progenitors with mpG/mBH > 
'nipc/'nii^S ^ 6 — 30. So the massloss is even more severe than for forming WDs. The progenitors evolve 
very quickly off the main sequence, and some explode as Type II SN, which immediately power a very fast 
wind. To keep the wind from severly damage a galaxy (cf. eq. 14 and eq. 6) would imply that these remnents 
make up a fraction 

Ibh < Ins < (0.1 - 1)% (20) 

of the total mass of a galaxy. In comparison, the upper limit on stellar BHs from microlensing is Jbh ^ 
(30— 100)% (Lasserre et al. 2000, Alcock et al. 2001). So our dynamical limit is much tighter. Some models 
of galaxy formations suggest a significant amount of massive black holes in galaxy halos, with black hole 
mass niBH ~ 10^ — 10^ Mq (Lacey & Ostriker 1985, Madau & Rees 2001). These black holes, however, must 
be pregalactic, in which case our dynamical limits do not apply. 

In summary, we show that intergalactic stars and globular clusters and the present day density profiles 
of galaxies offer a new diagnosis of early star bursts of massive stars in halos of galaxies, and the amount of 
remnents, such as halo WDs and BHs. A mild wind-induced expansion of the dark halos might have played 
some role in lowering the dense cold dark matter central density in dwarf galaxies (Gnedin & Zhao 2001) 
and high surface brightness galaxies as well. 

HSZ thanks Bernard Carr, Rodrigo Ibata, Mike Irwin, Gerry Gilmore, Donald Lynden-Bell, Jerry Ostriker, 
Tom Theuns and the anonymous referee for helpful comments. 



A. Appendix 

A.l. Effects of massloss on the size of the galaxy 

A rapid change of the potential of a galaxy after massloss can cause the galaxy to expand. The final size 
can be estimated with the virial theorem. Let Mj be the initial total mass, Mj = (1 — e)Mj be the final 
mass of the galaxy, and rj,e and r/,e be the initial and final effective sizes of the galaxy such that 

GM2 GM2 GM2 

^/^-V' (^1) 

where Wj, Wa, and W/ be the total potential energies of the galaxy before, right after, and long after the 
expansion. Then the total energies of the system before and after the expansion are given according to the 
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virial theorem by 



Wi + ^M.KL = \Wi, + ^M/y,L = \W}, (A2) 



where yj.jjjg is the root mean squared velocity of the system before or right after the instantaneous expansion. 
Combining the equations we find the overall expansion factor (cf. Hills 1980) 



fast 1 - 2e 



(A3) 



In comparision, if the niassloss is adiabatic, then a shell initially at radius maps to another shell at 
radius r adiabaticly conserving the angular momentum ^ G{Mq + Myj)ri = ^GM^r, so the expansion factor 



slow 



— • (A4) 



An interpolation of the equations for the fast and slow massloss can be made to model massloss of arbitrary 
speed. Generally the expansion factor X can be modeled as a function of the massloss fraction e as following 

X^!| = l + ^!2^, (A5) 

rf Coo - e 

where the constant Coo is the critical amount of massloss for the system to expands to infinity X ^ oo. For 
a wind of a range of speed we have 

[fast wind] 0.5 < Cqo < 1 [slow wind]. (A6) 

A. 2. Effects of massloss on the size of tracer orbits 

A rapid change of the potential of a galaxy after massloss can also alter the orbital distributions of its 
tracers, e.g., the halo globular clusters. The effects on the apocenter are easiest to work out for the tracers 
on initially circular orbits. A circular orbit of initial radius Oj can become radial because of the reduced 
restoring force after the massloss and reach a new apocenter r = a while keeping the pericenter at r = Oj. 
The new apocenter is determined by the specific energy equation, 



GM{r') 

2^ 



E = ^-1 '^dr (A7) 



where the first terms are the azimuthal kinetic energy for an orbit with the initial specific angular momentum 
Ji = GMi{ai)ai, the second terms on the r.h.s. of the equations are final potential of the galaxy, determined 
by the final distribution of the mass 

M(r) = (1 - e)Mi(rX-i) (A9) 

in the galaxy, where Mi(ri) is the initial mass distribution inside a radius r^, and X is the overall expansion 
factor. Approximate the initial potential by a scale- free potential with the mass profile given by 

Mi{n) = 5irf , > /3 < 3, (AlO) 

where /3 is the powerlaw slope, and 5i is a normalisation constant. Combining the above equations, we have 
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which relates between the expansion factor of the apocenter A — a/ai and the massloss fraction e in a fast 
wind; the factors A, X and e are the same everywhere for the assumed scale-free potential. For a Keplerian 
potential with (3 = 0, the above equation has a solution 



fast 1 - 2e 



^ (A12) 



It turns out the above solution is accurate for 10% for models with P ^ 0, and < e < 0.25. This applies 
to non-scalc-frcc potentials such as the Navarro, Prenk & White (1996) models as well because the result is 

insensitive to the mass profile parameter /?. 

In comparison, if the wind is slow, then the orbit will remain circular and expand adiabaticly by the same 
factor X as the galaxy, so 



= X = -^. (A13) 

slow ^ 



For wind of any speed, we can write 



(A14) 



Note A ^ X ^ 2 if lossing 50% of the orginal mass in a slow wind (eoo = 1), and ^ — > X — > oo if lossing 
the same amount of mass in a fast wind (coo = 0.5). 

A. 3. Effects of massloss on escaping stars 

A sudden weakening of the potential will boost the energy of some particles above their binding energy, 
and escape. The average specific energy of the system before and long-after the massloss is related to the 
characteristic potential energy (cf. eq.Al) by the virial theorem, so 

2Mi 4:ri,e' ^ 2Mf 4rf,e ^ ' 

Substitute in the expansion factor, we have 

Particles in a King model or any smoothly truncated model should be bound between two energies, 4>min at 
the bottom of the potential well, and (pmaK at the outer truncation. Long after the expansion, the truncation 
energy </)max is increased to 

= (1 - 2e)<^rr.ax. (A17) 



'^max 



During the sudden expansion, a particle at the bottom of the potential well will gain a potential energy 
— e^min without changing kinetic energy, hence the particle's energy E is increased to 

E' = E- ecl>^i^. (A18) 

So this particle could escape the final potential well if 

E' > 0;,,,. (A19) 
This suggests that all particles with initial energy E satisfying 

E - (pmaK > e(0min - St^^max) (A20) 



11 



could escape the final potential, i.e., the escaping particles occupy a narrow band at the upper end of the 
energy spectrum of the initial system, and the band width is oc e. The high energy tail of initial particles 
should follow a distribution function f{E) with 

^^^^^ oc f{E) oc ((/.^ax - E) (A21) 
to have a smooth truncation. So the amount of mass in the high energy tail M(> E) is given approximately 

by 

M(> E) (.^^ax - E)^ 



M{> Emin) (^i-max " ^min)^ ' 

Hence the escape fraction is given by 



(A22) 



g _ -^*,esc ^ ^2 ('/'min 2(^inax)'^ (A23) 
M*,esc + (?!>max - 0min)^ 

For a highly concentrated system, e.g., a King model with a core radius 100 times smaller than the truncation 
radius, we have 

l^maxl < |'/>min|- (A24) 



So for highly concentrated systems we have 

e* « e 



2 (A25) 



The above estimations can be generalized to winds of any speed with the following approximation 



^* ,esc 



M.esc + Mh. 



(1 - eoo)e 



(1 



2 



(A26) 



This parametrization is justified because it reproduces nicely two simple limiting cases: (i) no stars escape 
in a slow wind with Cqo ^ 1; (ii) all stars escape in a wind of maximum strength e — > eo 
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Gas outflow fraction e 




Gas outflow fraction e 

14 

Fig. 1. — Upper panel shows the halo expansion factor X — r jvi in case of a slow (dashed) or a fast wind 
(solid) lossing e fraction of the original mass of the galaxy everywhere. Middle panel shows the orbital 
expansion factor A = a /at vs. e. Bottom panel shows the fraction e, of all stars escaping from the galaxy. 




Fig. 2. — shows the distribution of the orbital parameters of halo globular clusters. Disk globulars are 
marked with starry symbols. The solid line indicates the predicted median apo-to-pericenter ratio in an 
isotropic model (van den Bosch et al. 1999), and the objects to the lower right of the dashed line have an 
extreme apo-to-peri ratio of > 10. 
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Fig. 3. — shows the effect of wind in the parameter space of MACHO-to-WIMP mass ratio fwo/i^ — fwo) 
vs. the progenitor-to-WD mass ratio mpc/mwD- Models between the two sohd hnes can have their halo size 
more than doubled by a fast wind. Models below the dashed line are allowed. The arrows show observational 
limits; the lower limits come from the number density by OHDHS (Oppenheimer et al. 2001), and IIBSG 
(Ibata et al. 2000); the upper limits come from the mass density by MACHO (Alcock et al. 1997, 2000) and 
EROS (Lasserre et al. 2000). A conversion scale between mpo and mwD is shown together with a cosmic 
"clock" marker indicating a 4000 Kelvin WD remnaj^ from a 10-Gyr-old burst; hotter WDs to the left of 
the marker, and cooler WDs to the right. 
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Fig. 4. — shows the temperature of WDs, which are predicted with the coohng isochrones of Richer et 
al. (2000) for hydrogen-rich WDs, together with a scale showing the initial-final mass conversion. Each 
isochronic curve is labeled with the age of the burst and corresponding redshift in bracket. 
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Fig. 5. — shows the aUowed halo WD and WIMP mix ratio fwo/i^ — fwo) and the temperature of WDs. 
The region above the sohd line is excluded because of too much mass loss from the WD progenitors which 
could double the size of the halo. The dashed line corresponds to a mild expansion of the halo. Here we take 
a 10 Gyr cooling isochronc and the initial- final relation from Richer et al. (2000, and references therein). The 
dynamical limit comes from eq. (14) and eq. (6). Several estimations for fwD from different observations 
are also roughly indicated by the labels and arrows at 4000 K, the typical temperature of observed WDs. 
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